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Abstract We established a method on measuring the D° - D° mixing parameter y for BESIII experiment 
at the BEPCII e + e~ collider. In this method, the doubly tagged ip(377Q) —* D°D° events, with one D decays 
to CP-eigenstates and the other D decays semileptonically, are used to reconstruct the signals. Since this 
analysis requires good e/n separation, a likelihood approach, which combines the dE/dx, time of flight and the 
electromagnetic shower detectors information, is used for particle identification. We estimate the sensitivity of 
the measurement of y to be 0.007 based on a 20/6 -1 fully simulated MC sample. 
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1 Introduction 



The mixing between a particle and its antiparticle 
has been observed experimentally in neutral K, B d 
and B a system. In the Standard Model, however, the 
mixing rate of the neutral D system is expected in 
general to be small and long-distance contributions 



make it difficult to be calculated^ 1 '— '— '^. Recently, 
several measurements^ — '— '— '^ present evidences for 
D° — D° mixing with the significance ranging from 
3 to 4 standard deviations. This highlights the need 
for independent measurements of the mixing parame- 
ters. Here, we present a study on measuring the mix- 
ing parameter y at BESIII experiment, which takes 
advantage of the correlated threshold production of 
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D° — D° pairs in e + e~ collisions. 

For the neutral D meson system, two mass eigen- 
states and flavor eigenstates are not equivalent and 
can be expressed as the following form of the two 
quantum states: 



\D A , B )=p\D°}±q\D°), 



(1) 



with eigenvalues of masses and widths to be m A}B and 
Y AB . Conventionally, the D°—D° mixing is described 
by two small dimcnsionless parameters: 



(2) 



Am _ Ar 

where Am = m A — m B , AT = T A — T B and T = 
(T A + T B )/2. The mixing rate R M is approximately 



R* 



■y 



(3) 



In the limit of CP conservation, the \D A ) and \D B ) 
denote the CP eigenstates. 

The mixing parameters can be measured in several 
ways. The i3-factories measured R M with semilep- 
tonic D° decay samples^ ^ . Reference^' ^ also 
gave an estimation on the sensitivities of R M mea- 
surement at BESIII. Other attempts'^'— '—'—'^ are 
based on the proper-time measurements of the neu- 
tral D meason decays. However, the time-dependent 
analyse are not possible at symmetric charm factory, 
which operates at the ^(3770) resonance. In this 
analysis, we utilize the quantum-coherent threshold 
production of D° — D° pairs in a state of definite 
C = — 1. Applying the kinematics of the process of 
e + e~ — > V(3770) — > D°D°, we can reconstruct both 
neutral D mesons (double tagging (DT) technique) to 
obtain clean samples to measure the mixing param- 
eters, the strong phase difference and the CP viola- 
tion. For the single D° meson decays into a CP eigen- 
state, the time-integrated decay rate can be written 



as 



[14. iS. 



To measure y at BESIII, only the electron chan- 
nels are used to reconstruct the semileptonic D° de- 
cays. In the muon channels, the transverse momen- 
tum of muon is too low to be efficiently identified by 
the BESIII muon detector. Thus, the e/ir separation 
plays an essential role to suppress the backgrounds. 
Fig. [U shows the momentum distribution of the elec- 
trons from the semileptonic D decays. The momen- 
tum distribution of the pions from s quark decays is 
similar to Fig. [T] As a result, the performance of elec- 
tron identification (e-ID) will determine the precision 
of the measurement of y parameter. 
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semileptonic decays 

The designed peak luminosity of BEPCII (Beijing 
Electron Position Collider) is 10 33 cm" 2 s _1 at beam 
energy E beam = 1.89 GeV, which is the highest in the 
tau-charm region ever planned and an unprecedented 
large number of ^(3770) events is expected. 

This paper is organized as follows: an improved 
electron identification technique for BESIII is de- 
scribed in Section 2. In Section 3, we describe the 
method on reconstructing the signals with Mente 
Carlo(MC) simulation samples. Section 4 presents 
the estimated sensitivity of y measurement. The sum- 
mary is presented in Section 5. 



T(D°^f CP± ) = 2A 2 CP± [lTy], (4) 2 Electron identification 



where fcp± is a CP eigenstate with eigenvalue ±1, 
and A C p± = \(fcp±\T~L\D°)\ is the magnitude of de- 
cay amplitude. If we consider the coherent _D-pair de- 
cays, in which one D decays into CP eigenstates and 
the other D decays semileptonically, the decay rate of 



(p°D° 



C=-l 



(I X) (fcp±) is described as 



ri6. 17. 181 



r l;C P = T[(l±X)(fcp)}nA% x A 2 CP , (5) 



where A t ± x = | {i ± -X"|7^|Z?°) |. Here, we neglect terms 
to order y 2 or higher since y is much smaller than 
unit. We, thus, can derive: 



1 



^l;CP+^CP~ 
^1;CP-LCP+ 



l;CP- 



^l;CP+^CP- 



(6) 



The BESIII detector operates at BEPCII and con- 
sists of a beryllium beam pipe, a helium-based small- 
celled drift chamber, Time-Of-Flight (TOF) coun- 
ters for particle identification, a CsI(Tl) crystal elec- 
tromagnetic calorimeter (EMC), a super-conducting 
solenoidal magnet with the field of 1 Tesla, and a 
muon identifier of Resistive Plate Counters (RPC) 
interleaved with the magnet yoke plates. The BESIII 
Offiline Software System (BOSS)^ of version 6.1.0 
is used for this analysis. The detector simulation' 2 ^ 
is based on GEANT4^. 

The BESIII detector has four subsystems for par- 
ticle identification: the dE/dx of the main drift 
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chamber(MDC), TOF, EMC and the muon counter. 
Among them, the dE/dx and the TOF systems are 
mainly used for hadron separation, the EMC pro- 
vides information for electron and photon identifi- 
cation, the MUC has good performance on muon 
identification^ . 

For electron identification, Refs^'^ illustrate 
the use of dE/dx of MDC and TOF information. 
Here, an improved e/n separation technique is intro- 
duced in the following sections. 

2.1 Electromagnetic calorimeter 

The BESIII electromagnetic calorimeter^^ is 
composed of one barrel and two endcap sections, cov- 
ering 93% of 4ir. There are a total of 44 rings of 
crystals along the z direction in the barrel, each with 
120 crystals. And there are 6 layers in the endcap, 
with different number of crystals in each layer. The 
entire calorimeter has 6240 CsI(Tl) crystals with a 
total weight of about 24 tons. The energy resolution 
is expected to be 2.5% and the spatial resolution is 
expected to be 0.6 cm for 1 GeV/c photon. 

The primary function of the EMC is to precisely 
measure the energies and positions of electron and 
photon. In order to distinguish electron from hadron, 
we make use of significant differences in energy depo- 
sition and the shower shape of different type of the 
particles. 

2.2 Variables used in e-ID 

The following variables are used to identify the 
electron from pion: 

1) Ratio of the energy measured by the EMC and the 

momentum of the charged track by the MDC 
(E/P). 

Ratio of the energy measured by the EMC and the 
momentum of the charged track by the MDC (E/p). 
When an electron passes through the calorimeter, the 
electron produces electromagnetic shower and loses 
its energy by pair-production, Bremsstrahlung and 
ionizing/exciting atomic electrons. Since the mass of 
electron is negligible in the energy range of interest, 
we expect to have the ratio E/p=l within the mea- 
surement errors. For hadrons, the E/p is typically 
smaller than one. 

2) Lateral shower shape at the EMC. 

In order to enhance the separation between the 
electrons and the interacting hadrons, the lateral 
shower shape can also be utilized. These vari- 
ables include: E aeed /E 3x3 ,E 3x3 /E 5x5 and the second- 
moment. Here the E seed is the energy deposited in 
the central crystal, the -£3x3 and E 5x5 represent the 



energy deposit in the 3x3 and 5x5 crystal array, 
respectively. The second-moment S is defined as 



S = 



(7) 



where Ei is the energy deposit in the i-th crystal, 
and di is the distance between the i-th crystal and 
the center position of reconstructed shower. Detailed 
description of E/p and the lateral shower shape can 
be found in Ref ^ . 

3) Longitudinal shower shape at the EMC. 

The longitudinal shower shape provides additional in- 
formation for electron identification. The variable 
A<f>, between the polar angles where the track inter- 
sects the EMC and the shower center, can be used. 
The distributions of A(f> for electron and pion are 
drawn in Fig. [2l The center of electron showers is 
closer to the impact point of track on EMC since the 
electron showers reach their maximum earlier than 
hadrons. 





Fig. 2. A<fi of (a) electron (b) pion. 
2.3 The correlation between variables 

The E/p ratio, lateral shower shape and longitu- 
dinal shower shape are all depending on the deposited 
energy in the crystals. Thus, these variables may be 
correlated. We calculate the correlation coefficients 
Pij between the E/p, E 3x3 / E 5x5 and A(f> using the 
function: 



i,3 



- Xj ) , Pij 



M, 



y/Mu X Mjj ' 

(8) 

where i,j are the indices of the variable names. Fig- 
ure. [3] shows the correlation between any two of the 
variables of electron and pion, respectively, with the 
momentum ranging from 0.2 GeV/c to 2.0 GeV/c. 
Here, the x-axis represents the particle momentum 
and the y-axis represents the correlation coefficient 
Pij. The distribution indicates strong correlation be- 
tween the variables. 
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Fig. 3. Correlations between (a)E/p and 
E aee d / E-ixz of electron; (b)E/p and 
Esx3/E 5 yj, of electron; {c)E seed / E Sx :>, 
and £3x3/^5x5 of electron; (d)E/p and 
E seed /E3x3 of pion; (e)E/p and E 3x3 /E 5x5 of 
pion; (f)-E see d/i?3x3 and E 3x3 /E 5x5 of pion. 



Fig. 4. The NN outputs of (a) pion (b) electron 
samples with the momentum ranging from 
0.2GeV/c to 1.6GeV/c 

2.5 The performance check 

To combine the dE/dx, TOF and EMC informa- 
tion, the likelihood approach^ is adopted. Firstly, 
the likelihood value of each subsystem is calculated. 
Then, the total likelihood value of each hypothesis is 
calculated by the following formula: 



2.4 PID Algorithm 

Considering the correlations between the vari- 
ables, the traditional method for particle identifica- 
tion may be underperforming. In the e-ID, we im- 
plement the artificial neural network (ANN) «a to 
provide a general framework for estimating non-linear 
functional mapping between the input variables and 
the output variable. For the neural network (NN) 
training, we use the momentum, traverse momentum 
and other six discriminants (total deposit energy Es- 
eed, E3x3, E5x5, second moment and A<f>) as the in- 
put variables. The network we choose has one hidden 
layer with 16 neurons and one output value. Fig- 
ure. 0] shows the two-dimension distributions of the 
output value versus the momentum of the electrons 
and pions. It is obvious that the distribution of the 
output value depends on the momentum, especially 
at low momentum region. Thus, it is unsuitable to 
apply a single cut on the output value to separate the 
electrons from the pions. In practice, we construct 
probability density function (PDF) of the NN output 
value at every 0.1 GeV/c momentum bin. The PDF 
is obtained from fitting the nearest 4 bins of the NN 
output value, with the third-order polynomial func- 
tion. Then, the PDF value of the NN output can be 
extracted from the fit. Finally, we make the PID de- 
cision by comparing the likelihood values of electron 
and pion hypothesis. 



L to t — L 



dE/dx 



(9) 



where L 



dE/dx 



and L TOF represent the likelihood value 
of dE/dx and TOF subsystems respectively. Finally 
the likelihood ratio of electron hypothesis is defined 
as: 

lhf e = -r^i-, (10) 



where L e and L n are the total likelihood value of elec- 
tron and 7r hypothesis. To check the performance of 
the e/ir separation, both the electron and pion sam- 
ples are generated with the momentum ranging from 
0.2 Gev/c to 1.6 GeV/c, by using single particle gen- 
erator. Fig. EKa) shows the electron likelihood ra- 
tio distributions of these samples. For a particle to 
be identified as an electron, we require lhf e > 0.5. 
Fig. 0(b) shows the combined e /w separation perfor- 
mance using the dE/dx, TOF and EMC systems. 
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3 Simulation and reconstruction 

3.1 The reconstruction of CP tags 

For the neutral D meson decays, the main decay 
modes of CP+ eigenstate are K + K~ , tt + it~ , K s ir°ir°, 
7T 7r°, K S K S and p°7r°. The CP— eigenstates decay 
through the modes K s (7r° , p° ,r),r)' ,</>>&)• Considering 
the branching ratio and the reconstruction efficiency, 
we only simulated the K + K~, ir + ir~ for CP+ tag- 
ging, and the K s (n° 7 i],f]') for CP— tagging. 

For selecting the charged tracks, the following se- 
lection criteria are adopted: 

1) All charged tracks must have a good helix fit, and 

are required to be measured in the fiducial re- 
gion of MDC; 

2) Their parameters must be corrected for energy 

loss and multiple scattering according to the as- 
signed mass hypotheses; 

3) The tracks not associated with K° s reconstruction 

are required to be originated from the interac- 
tion point(IP). 

For reconstructing the CP+ eigenstates, two 
opposite-charged tracks of K or 7r are selected with 
the requirements that they are from IP and to pass 
a common vertex constraint. To identify a track as a 
7r or K, we use the likelihood method to combine the 
information of dE/dx and TOF with the likelihood 
fraction of 7r or K greater than 0.5. Then the beam 
constrained mass(Mi, c ) of the D meson is used to dis- 
tinguish the signal and background, and it is defined 
as: 

M bc EE ^El am - (]>>) 2 = \fKZn - {Pd)\ (11) 

where the E beam is the beam energy, the Pi is the mo- 
mentum of the i-th track and the p D = ^2,Pi is the 
momentum of the reconstructed D meson. 

For tagging the CP— eigenstates, we need to re- 
construct the neutral mesons K Sl 7r°, 77 and 77 . The 
K s candidates are reconstructed through the decay 
of Ks — > 7r + 7r~. The decay vertex formed by 7r + 7r~ 
pair is required to be away from the interaction point, 
and the momentum vector of 7r + 7r~ pair must be 
aligned with the position vector of the decay ver- 
tex to the IP. Here we set L vtx /a vtx to be greater 
than 2, where L vtx and a vtx are the measured decay 
length and the error of the decay length of the K s . 
The tt + tt~ invariant mass is required to be consis- 
tent with the Kg nominal mass within ±10 MeV. 
To identify the neutral tracks, one has to address 
a number of processes which can produce both real 
and spurious showers in EMC. The major source of 



these "fake photons" arises from hadronic interaction, 
which can create a "split-off" shower. This shower 
does not associate with the main shower and may be 
recognized as a photon. Other sources of fake pho- 
tons include particle decays, back splash, beam as- 
sociated background and electronic noise. To reject 
"fake photons", the selection criteria for "good pho- 
ton" include a deposit energy cut, and a spatial cut, 
which requires that the cluster is isolated from the 
nearest charged tracks. These "cuts" are set to be 
E 1 > 40MeV and A C7 > 18°, where E 1 and A C7 rep- 
resent the deposited energy and the crossing angle of 
the cluster to the nearest charged track, respectively. 
The neutral pions are reconstructed from 7r° — > 77 
decays using the photons observed in the barrel and 
endcap regions of EMC. At the energies of interest, 
a 7T° decays into two isolated photons. In addition, 
we also reconstruct 77/77 candidates in the modes of 
77 — > 77,77 — > 7r + 7r~7r°,?7 — > 7/5 and 77 — >?77r + 7r~. For 
these modes, 3<r consistency with the ir°/r]/ri mass 
is required, followed by a kinematic mass constraint. 
For CP— eigenstates, the beam constrained mass is 
also used to select the signal. 

Under the environment of BOSS 6.1.0, we simu- 
lated D° — D° pairs production at the -0(3770) peak 
with one D decayed into CP eigenstates and the other 
D decayed semileptonically. The CP+ eigenstates 
are decayed through tt + it~ and K + K~ according to 
their branching ratios. For CP— eigenstates, the de- 
cay modes K s ir°, K s r\ and K s r\ are included. In the 
K Sl 7T°, 77 and 77' decays, the decay modes are listed as 

follows: Kg — * 7T + 7T~ , 7T° — > 77, T] — > 77, 7] —> 7T + 7T _ 7T° , 

77 ^7p°, and 77 — >r)Tr + n~. For the CP+ and CP— 
eigenstates, we generated 30,000 events for each MC 
sample. The distributions of the beam constrained 
mass of D meson are shown in Fig. [6] 
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Fig. 6. The M bc of (a)the CP- tags; (b)the 
CP+ tags. 

3.2 The reconstruction of semileptonic tags 

For tagging the semileptonic decays, we use the 
decay mode D° — > K~e + v e . To reconstruct the neu- 
tral D meson, good tracks for one electron and one 
kaon candidate are required. The good track selec- 
tion criteria are the same as the CP tagging discussed 
in Section 3.1. The electron and kaon candidates are 
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also required to be from the IP, and the likelihood ra- 
tio of the electron and kaon must both be greater than 
0.5. Moreover, the two tracks need to pass a common 
vertex constraint. After the electron and pion selec- 
tions, a standard partial reconstruction technique is 
applied to this semileptonic decay channel with one 
neutrino associated. Here, we use the "missing mass" 
{Umiss) of neutrino to select the signal candidates. 
The "missing mass" is defined as follows: 



■ E„ 



Prr 



(12) 



where E misa = E Dtag - E K - E e and P miss = \p Dtag - 
p K — p e | are the missing energy and momentum of 
the neutrino. Here, E K e and p K e are the measured 
energy and momentum of the selected kaon and elec- 
tron track. E Dtag is the energy of the D meson, which 
is equal to the beam energy. Po tag — —Pd cp i s the 
3-momentum of the D meson, which can be obtained 
from the reconstructed momentum of the CP tagged 
D meson. For neutrino, the energy and momentum 
are equal. Thus, the distribution of U misa must have 
a mean value at zero. The distribution of U misa is 
shown in Fig. [71 We apply a 3er cut on the U mias to 
select the semileptonic D decays. 
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Fig. 7. Umiss of v e for (a) CP+ tags and (b) 
CP- tags. 



4 The sensitivity of y 

Table Q] shows the reconstruction efficiency and 
the number of estimated doubly-tagged events for dif- 
ferent simulated decay channels. For ~ 20fb _1 lumi- 
nosity at -0(3770) peak, which approximately corre- 
sponds to four years data taking at BESIII, about 
8.0 x 10 7 D° — D° pairs can be produced. Accord- 
ing to the full simulation, about 11000 doubly tagged 
CP+ decays and 9000 doubly tagged CP— decays 
can be reconstructed. 



Table 1. The efficiency and expected events for 
8.0 x 10 7 D° - D° decays for different decay 
channels. 



decay mode 


efficiency 


event estimation 


K~e + v e 
K-K+,-n-n+ 


40% 


11701 


K~e + v e 
K s n° 


16.8% 


7345 


K~e + v e 
KsV 


7.7% 


715 


K~e + v e 
K sV ' 


4.7% 


953 



For a small y, to calculate the a y of Equation(6), we 
ignore the statistical error from single tagged events. 
Hence, the statistical error of y parameter can be ob- 
tained from the following equation: 



1 



1 1 

N~i + AV 



(13) 



where N 1 and N 2 represent the reconstructed doubly- 
tagged CP+ and CP— events. As a result, the a v is 
estimated to be 0.007 with - 20fb _1 data at 0(3770) 
peak in this analysis. Since the double tagging tech- 
nique is adopted here, the background effect can be 
ignored comparing to the statistical sensitivity esti- 
mated above. 



5 Summary 

In this paper, we presented a MC study on mea- 
suring the D°—D° mixing parameter y at the BESIII 
experiment. Based on a 20/6~ 1 fully simulated MC 
sample of -0(3770) resonance decays, we estimated 
the sensitivity of the y measurement to be 0.007. In 
this analysis, the double tagging technique was used 
for reconstructing the D meson pairs. Here, the sig- 
nal is reconstructed such that one D decays to CP 
eigenstates and the other D decays semileptonically. 
The electron identification is essential for this analy- 
sis. We improved the e-ID technique for BESIII ex- 
periment, which can also be applied to many other 
important physics topics. Our next step is to include 
more semileptonic decay modes, such as D° — > K*ev R , 
into this analysis to improve the sensitivity of y mea- 
surement. 
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Abstract We established a method on measuring the D° — D° mixing parameter y for BESIII experiment 
at the BEPCII e + e~ collider. In this method, the doubly tagged -0(3770) -^D a D° events, with one D decays 
to CP-eigenstates and the other D decays semileptonically, are used to reconstruct the signals. Since this 
analysis requires good e/n separation, a likelihood approach, which combines the dE/dx, time of flight and the 
electromagnetic shower detectors information, is used for particle identification. We estimate the sensitivity of 
the measurement of y to be 0.007 based on a 20/6 -1 fully simulated MC sample. 
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1 Introduction 

The mixing between a particle and its antiparticle 
has been observed experimentally in neutral K, B d 
and B a system. In the Standard Model, however, the 
mixing rate of the neutral D system is expected in 
general to be small and long-distance contributions 
make it difficult to be calculated^'— '—'^ . Recently, 
several measurements^'—'—'—'^ present evidences for 
D° — D° mixing with the significance ranging from 
3 to 4 standard deviations. This highlights the need 
for independent measurements of the mixing parame- 
ters. Here, we present a study on measuring the mix- 
ing parameter y at BESIII experiment, which takes 



advantage of the correlated threshold production of 
D° — D° pairs in e + e~ collisions. 

For the neutral D meson system, two mass eigen- 
states and flavor eigenstates are not equivalent and 
can be expressed as the following form of the two 
quantum states: 



\D A>B )=p\D°)±q\D°), 



(1) 



with eigenvalues of masses and widths to be mj ]B and 
r.4,s- Conventionally, the D°—D° mixing is described 
by two small dimcnsionlcss parameters: 

_ Am _ AT 
x = —,y = —, 

where Am = m, A — m B , AT = T A — T B and T 



(2) 



Received 29 May 2008, Revised 11 June 2008 
* Supported by National Natural Science Foundation of China (10491300,10491303,10735080), Research and Development 
Project of Important Scientific Equipment of CAS (H7292330S7), 100 Talents Programme of CAS (U-25, U-54,U-612) and Scientific 
Research Fund of GUCAS(110200M202) 

1) E-mail: huangb@ihep.ac.cn 

2) E-mail: zhengyh@gucas.ac.cn 

3) E-mail: liwd@ihep.ac.cn 

i — m 



No. 



2 



(T A + T B )/2. The mixing rate R M is approximately 

x 2 + y 2 



500 



(3) 



In the limit of CP conservation, the \D A ) and \D B ) 
denote the CP eigenstates. 

The mixing parameters can be measured in several 
ways. The i?-factories measured R M with semilep- 
tonic D° decay samples^ Reference^' ^ also 
gave an estimation on the sensitivities of R M mea- 
surement at BESIII. Other attempts'^'— '—'—'^ are 
based on the proper-time measurements of the neu- 
tral D meason decays. However, the time-dependent 
analyse are not possible at symmetric charm factory, 
which operates at the i/"(3770) resonance. In this 
analysis, we utilize the quantum-coherent threshold 
production of D° — D° pairs in a state of definite 
C = — 1. Applying the kinematics of the process of 
e + e~ — > ^(3770) — > D°D°, we can reconstruct both 
neutral D mesons (double tagging (DT) technique) to 
obtain clean samples to measure the mixing param- 
eters, the strong phase difference and the CP viola- 
tion. For the single D° meson decays into a CP eigen- 
state, the time-integrated decay rate can be written 



as 



[14. iS. 



rv 



T(D ^f CP± ) = 2A 2 CP± [l T y}, (4) 



where fcp± is a CP eigenstate with eigenvalue ±1, 
and A C p± = \{fcp±\7~L\D°)\ is the magnitude of de- 
cay amplitude. If we consider the coherent D-pair de- 
cays, in which one D decays into CP eigenstates and 
the other D decays semileptonically, the decay rate of 
(D°D^) C ^ 1 -» (/±X) (/ OF± ) is described as^'^'^: 



T l , CP = T[{l ± X){f GP )]; 



• A 2 , A 2 



(5) 



where A t ± x = K^^T |W|£) )|. Here, we neglect terms 
to order y 2 or higher since y is much smaller than 
unit. We, thus, can derive: 



i /r 



y = 



l;CP+ L CP- 



l;CP- L CF^ 



4 vr 



r 



l;CP+ L CP- 



(6) 



To measure y at BESIII, only the electron chan- 
nels are used to reconstruct the semileptonic D° de- 
cays. In the muon channels, the transverse momen- 
tum of muon is too low to be efficiently identified by 
the BESIII muon detector. Thus, the e/n separation 
plays an essential role to suppress the backgrounds. 
Fig. [1] shows the momentum distribution of the elec- 
trons from the semileptonic D decays. The momen- 
tum distribution of the pions from s quark decays is 
similar to Fig. [T] As a result, the performance of elec- 
tron identification (c-ID) will determine the precision 
of the measurement of y parameter. 




0.4 0.6 0.8 

P/(GeV/c) 

Fig. 1. Momentum of the electron from D° 
semileptonic decays 



The designed peak luminosity of BEPCII (Beijing 
Electron Position Collider) is 10 33 cm" 2 s _1 at beam 
energy E beam = 1.89 GeV, which is the highest in the 
tau-charm region ever planned and an unprecedented 
large number of ifj{?>77Q) events is expected. 

This paper is organized as follows: an improved 
electron identification technique for BESIII is de- 
scribed in Section 2. In Section 3, we describe the 
method on reconstructing the signals with Mente 
Carlo(MC) simulation samples. Section 4 presents 
the estimated sensitivity of y measurement. The sum- 
mary is presented in Section 5. 

2 Electron identification 

The BESIII detector operates at BEPCII and con- 
sists of a beryllium beam pipe, a helium-based small- 
celled drift chamber, Time-Of-Flight (TOF) coun- 
ters for particle identification, a CsI(Tl) crystal elec- 
tromagnetic calorimeter (EMC), a super-conducting 
solenoidal magnet with the field of 1 Tesla, and a 
muon identifier of Resistive Plate Counters (RPC) 
interleaved with the magnet yoke plates. The BESIII 
Offiline Software System (BOSS)^ of version 6.1.0 
is used for this analysis. The detector simulation^ 
is based on GEANT4^ . 

The BESIII detector has four subsystems for par- 
ticle identification: the dE/dx of the main drift 
chamber (MDC), TOF, EMC and the muon counter. 
Among them, the dE/dx and the TOF systems are 
mainly used for hadron separation, the EMC pro- 
vides information for electron and photon identifi- 
cation, the MUC has good performance on muon 
identification^ . 

For electron identification, Refs^'^ illustrate 
the use of dE/dx of MDC and TOF information. 
Here, an improved e/ir separation technique is intro- 
duced in the following sections. 
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2.1 Electromagnetic calorimeter 

The BESIII electromagnetic calorimeter^^ is 
composed of one barrel and two endcap sections, cov- 
ering 93% of 4ir. There are a total of 44 rings of 
crystals along the z direction in the barrel, each with 
120 crystals. And there are 6 layers in the endcap, 
with different number of crystals in each layer. The 
entire calorimeter has 6240 CsI(Tl) crystals with a 
total weight of about 24 tons. The energy resolution 
is expected to be 2.5% and the spatial resolution is 
expected to be 0.6 cm for 1 GeV/c photon. 

The primary function of the EMC is to precisely 
measure the energies and positions of electron and 
photon. In order to distinguish electron from hadron, 
we make use of significant differences in energy depo- 
sition and the shower shape of different type of the 
particles. 

2.2 Variables used in e-ID 

The following variables are used to identify the 
electron from pion: 

1) Ratio of the energy measured by the EMC and the 

momentum of the charged track by the MDC 
(E/p). 

Ratio of the energy measured by the EMC and the 
momentum of the charged track by the MDC (E/p). 
When an electron passes through the calorimeter, the 
electron produces electromagnetic shower and loses 
its energy by pair-production, Bremsstrahlung and 
ionizing/exciting atomic electrons. Since the mass of 
electron is negligible in the energy range of interest, 
we expect to have the ratio E/p=l within the mea- 
surement errors. For hadrons, the E/p is typically 
smaller than one. 

2) Lateral shower shape at the EMC. 

In order to enhance the separation between the 
electrons and the interacting hadrons, the lateral 
shower shape can also be utilized. These vari- 
ables include: E st , ed /E 3x 3 , E 3x3 / E 5x5 and the second- 
moment. Here the E aeed is the energy deposited in 
the central crystal, the E 3x3 and E 5xS represent the 
energy deposit in the 3x3 and 5x5 crystal array, 
respectively. The second-moment S is defined as 



S = 



(7) 



where E t is the energy deposit in the i-th crystal, 
and di is the distance between the i-th crystal and 
the center position of reconstructed shower. Detailed 
description of E/p and the lateral shower shape can 
be found in Ref 23. 



3) Longitudinal shower shape at the EMC. 

The longitudinal shower shape provides additional in- 
formation for electron identification. The variable 
A(j), between the polar angles where the track inter- 
sects the EMC and the shower center, can be used. 
The distributions of A<f> for electron and pion are 
drawn in Fig. [2l The center of electron showers is 
closer to the impact point of track on EMC since the 
electron showers reach their maximum earlier than 
hadrons. 





Fig. 2. A<f> of (a) electron (b) pion. 

2.3 The correlation between variables 

The E/p ratio, lateral shower shape and longitu- 
dinal shower shape are all depending on the deposited 
energy in the crystals. Thus, these variables may be 
correlated. We calculate the correlation coefficients 
Pij between the E/p, E 3x3 / E 5x5 and A(f> using the 
function: 



Mi. 



y/M lt X Mjj ' 

(8) 

where i, j are the indices of the variable names. Fig- 
ure. [3] shows the correlation between any two of the 
variables of electron and pion, respectively, with the 
momentum ranging from 0.2 GeV/c to 2.0 GeV/c. 
Here, the x-axis represents the particle momentum 
and the y-axis represents the correlation coefficient 
p^. The distribution indicates strong correlation be- 
tween the variables. 



(a) 



(c) 



(b) 
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Fig. 3. Correlations between (a)E/p and 
E seed /E 3x3 of electron; (b)E/p and 
E 3x3 /E 5x5 of electron; (c)E seed /E 3x3 
and E 3x3 /Es X 5 of electron; (d)E/p and 
E seed /E3x3 of pion; (e)E/p and E 3x3 /E 5x5 of 
pion; (i)E seed /E 3x3 and E 3x3 /E 5x5 of pion. 
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2.4 PID Algorithm 

Considering the correlations between the vari- 
ables, the traditional method for particle identifica- 
tion may be underperforming. In the e-ID, we im- 



[20] 



to 



plement the artificial neural network (ANN) 
provide a general framework for estimating non-linear 
functional mapping between the input variables and 
the output variable. For the neural network (NN) 
training, we use the momentum, traverse momentum 
and other six discriminants (total deposit energy, Es- 
eed, E3x3, E5x5, second moment and A<f>) as the in- 
put variables. The network we choose has one hidden 
layer with 16 neurons and one output value. Fig- 
ure. 0] shows the two-dimension distributions of the 
output value versus the momentum of the electrons 
and pions. It is obvious that the distribution of the 
output value depends on the momentum, especially 
at low momentum region. Thus, it is unsuitable to 
apply a single cut on the output value to separate the 
electrons from the pions. In practice, we construct 
probability density function (PDF) of the NN output 
value at every 0.1 GeV/c momentum bin. The PDF 
is obtained from fitting the nearest 4 bins of the NN 
output value, with the third-order polynomial func- 
tion. Then, the PDF value of the NN output can be 
extracted from the fit. Finally, we make the PID de- 
cision by comparing the likelihood values of electron 
and pion hypothesis. 




ANN output 

Fig. 4. The NN outputs of (a) pion (b) electron 
samples with the momentum ranging from 
0.2GeV/c to 1.6GeV/c 

2.5 The performance check 

To combine the dE/dx, TOF and EMC informa- 
tion, the likelihood approach^ is adopted. Firstly, 
the likelihood value of each subsystem is calculated. 
Then, the total likelihood value of each hypothesis is 
calculated by the following formula: 



where L dE / dx and L TOF represent the likelihood value 
of dE/dx and TOF subsystems respectively. Finally, 
the likelihood ratio of electron hypothesis is defined 
as: 

lhfe = -r^i-, (10) 



where L e and L n are the total likelihood value of elec- 
tron and 7r hypothesis. To check the performance of 
the e/n separation, both the electron and pion sam- 
ples are generated with the momentum ranging from 
0.2 Gev/c to 1.6 GeV/c, by using single particle gen- 
erator. Fig. |3fa) shows the electron likelihood ra- 
tio distributions of these samples. For a particle to 
be identified as an electron, we require lhf e > 0.5. 
Fig. Efb) shows the combined e /it separation perfor- 
mance using the dE/dx, TOF and EMC systems. 




(b) 



JT 



P /(Gev/c) 

Fig. 5. {a)lhf e of electron and n samples; (b) 
performance of e/n seperation. 



3 Simulation and reconstruction 



3.1 The reconstruction of CP tags 



Ltot — L. 



dE/dx * ^TOF 



(9) 



For the neutral D meson decays, the main decay 
modes of CP+ eigenstate are K + K~ , 7r + 7r~, i^ S 7r°7r°, 
7r°7r°, K S K S and p°tt°. The CP— eigenstates decay 
through the modes K s (ir° ,p° ,77,77' ,<fi,uj). Considering 
the branching ratio and the reconstruction efficiency, 
we only simulated the K + K~, ir + ir~ for CP+ tag- 
ging, and the K s (tt° ,rj,r}') for CP— tagging. 

For selecting the charged tracks, the following se- 
lection criteria are adopted: 
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1) All charged tracks must have a good helix fit, and 

are required to be measured in the fiducial re- 
gion of MDC; 

2) Their parameters must be corrected for energy 

loss and multiple scattering according to the as- 
signed mass hypotheses; 

3) The tracks not associated with Kg reconstruction 

are required to be originated from the interac- 
tion point(IP). 

For reconstructing the CP+ eigenstates, two 
opposite-charged tracks of K or 7r are selected with 
the requirements that they are from IP and to pass 
a common vertex constraint. To identify a track as a 
7r or K, we use the likelihood method to combine the 
information of dE/dx and TOF with the likelihood 
fraction of it or K greater than 0.5. Then the beam 
constrained mass(M i)c ) of the D meson is used to dis- 
tinguish the signal and background, and it is defined 
as: 

M bc = y/E* eam - (]>>) 2 = y/El am - {p D f, (11) 

where the E beam is the beam energy, the Pi is the mo- 
mentum of the i-th track and the p D = ^2,Pi is the 
momentum of the reconstructed D meson. 

For tagging the CP— eigenstates, we need to re- 
construct the neutral mesons K Sl 7r°, 77 and 77 . The 
K s candidates are reconstructed through the decay 
of K s — > tt + t:~ . The decay vertex formed by 7r + 7r _ 
pair is required to be away from the interaction point, 
and the momentum vector of tt + tt~ pair must be 
aligned with the position vector of the decay ver- 
tex to the IP. Here we set L vtx /a vtx to be greater 
than 2, where L vtx and a vtx are the measured decay 
length and the error of the decay length of the K s . 
The tt + tt~ invariant mass is required to be consis- 
tent with the Kg nominal mass within ±10 MeV. 
To identify the neutral tracks, one has to address 
a number of processes which can produce both real 
and spurious showers in EMC. The major source of 
these "fake photons" arises from hadronic interaction, 
which can create a "split-off" shower. This shower 
does not associate with the main shower and may be 
recognized as a photon. Other sources of fake pho- 
tons include particle decays, back splash, beam as- 
sociated background and electronic noise. To reject 
"fake photons" , the selection criteria for "good pho- 
ton" include a deposit energy cut, and a spatial cut, 
which requires that the cluster is isolated from the 
nearest charged tracks. These "cuts" are set to be 
E 1 > iOMeV and A C7 > 18°, where E 1 and A C7 rep- 
resent the deposited energy and the crossing angle of 
the cluster to the nearest charged track, respectively. 



The neutral pions are reconstructed from ir° — > 77 
decays using the photons observed in the barrel and 
endcap regions of EMC. At the energies of interest, 
a 7T° decays into two isolated photons. In addition, 
we also reconstruct 77/77 candidates in the modes of 
77 — > 77,77 — > 7r + 7r _ 7r°,77 —> 7/3 and 77 — >r)Tr + n~. For 
these modes, 3<r consistency with the ir°/r]/r] mass 
is required, followed by a kinematic mass constraint. 
For CP— eigenstates, the beam constrained mass is 
also used to select the signal. 

Under the environment of BOSS 6.1.0, we simu- 
lated D° — D° pairs production at the -0(3770) peak 
with one D decayed into CP eigenstates and the other 
D decayed semileptonically. The CP+ eigenstates 
are decayed through tt + it~ and K + K~ according to 
their branching ratios. For CP— eigenstates, the de- 
cay modes K s ir°, K s r\ and K s r\ are included. In the 
K s , 7T°, 77 and 77 decays, the decay modes are listed as 

follows: K S — > 7T + 7T~, 7T°^77, 77^77, 77 — > 7T + 7T~7r , 

77 — >"fp°, and 77 ^ 777r + 7r". For the CP+ and CP— 
eigenstates, we generated 30,000 events for each MC 
sample. The distributions of the beam constrained 
mass of D meson are shown in Fig. [6] 







\ (b) 









1.84 1.85 1.S6 1.87 1.S8 1.S9 1.85 1.855 1.86 1.865 1.87 1.875 1.88 

Mbc/(GeV/c I ) Mte/(GeV/c 2 ) 
Fig. 6. The M bc of (a)the CP- tags; (b)the 
CP+ tags. 

3.2 The reconstruction of semileptonic tags 

For tagging the semileptonic decays, we use the 
decay mode D° — > K~e + v e . To reconstruct the neu- 
tral D meson, good tracks for one electron and one 
kaon candidate are required. The good track selec- 
tion criteria are the same as the CP tagging discussed 
in Section 3.1. The electron and kaon candidates are 
also required to be from the IP, and the likelihood ra- 
tio of the electron and kaon must both be greater than 
0.5. Moreover, the two tracks need to pass a common 
vertex constraint. After the electron and pion selec- 
tions, a standard partial reconstruction technique is 
applied to this semileptonic decay channel with one 
neutrino associated. Here, we use the "missing mass" 
(Umiss) of neutrino to select the signal candidates. 
The "missing mass" is defined as follows: 

Urniss = E m ^ ss Pmiss ) (^^) 

where E miss = E Dtag - E K - E e and P misa = \p Dtag - 
p K — p e | are the missing energy and momentum of 
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the neutrino. Here, E K<£ and p x , e are the measured 
energy and momentum of the selected kaon and elec- 
tron track. Eo t is the energy of the D meson, which 
is equal to the beam energy. PD tag — ~Pd cp is the 
3-momentum of the D meson, which can be obtained 
from the reconstructed momentum of the CP tagged 
D meson. For neutrino, the energy and momentum 
are equal. Thus, the distribution of U mi3S must have 
a mean value at zero. The distribution of U misa is 
shown in Fig. [71 We apply a 3cr cut on the U miS3 to 
select the semileptonic D decays. 




Uiiss/GeV U.i„/GeV 



Fig. 7. Umisa of v e for (a) CP+ tags and (b) 
CP- tags. 
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For a small y, to calculate the a y of Equation(6), we 
ignore the statistical error from single tagged events. 
Hence, the statistical error of y parameter can be ob- 
tained from the following equation: 

(t v = -x x — + — , (13) 

where N 1 and N 2 represent the reconstructed doubly- 
tagged CP+ and CP— events. As a result, the a y is 
estimated to be 0.007 with ~ 20fb~ 1 data at V(3770) 
peak in this analysis. Since the double tagging tech- 
nique is adopted here, the background effect can be 
ignored comparing to the statistical sensitivity esti- 
mated above. 



4 The sensitivity of y 

Table Q] shows the reconstruction efficiency and 
the number of estimated doubly-tagged events for dif- 
ferent simulated decay channels. For ~ 20fb _1 lumi- 
nosity at i/>(3770) peak, which approximately corre- 
sponds to four years data taking at BESIII, about 
8.0 x 10 7 D° — D° pairs can be produced. Accord- 
ing to the full simulation, about 11000 doubly tagged 
CP+ decays and 9000 doubly tagged CP— decays 
can be reconstructed. 

Table 1. The efficiency and expected events for 
8.0 x 10 7 D° -D° decays for different decay 
channels. 



5 Summary 

In this paper, we presented a MC study on mea- 
suring the D°—D° mixing parameter y at the BESIII 
experiment. Based on a 20/6~ 1 fully simulated MC 
sample of i/>(3770) resonance decays, we estimated 
the sensitivity of the y measurement to be 0.007. In 
this analysis, the double tagging technique was used 
for reconstructing the D meson pairs. Here, the sig- 
nal is reconstructed such that one D decays to CP 
eigenstates and the other D decays semileptonically. 
The electron identification is essential for this analy- 
sis. We improved the e-ID technique for BESIII ex- 
periment, which can also be applied to many other 
important physics topics. Our next step is to include 
more semileptonic decay modes, such as D° — > K'ev,,, 
into this analysis to improve the sensitivity of y mea- 
surement. 
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